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Abstract 

We consider membranes interacting via short, intermediate and 
long stickers. The effects of the intermediate stickers on the lateral 
phase separation of the membranes are studied via mean-field approx- 
imation. The critical potential depth of the stickers increases in the 
presence of the intermediate sticker. The lateral phase separation of 
the membrane thus suppressed by the intermediate stickers. Consid- 
ering membranes interacting with short and long stickers, the effect 
of confinement on the phase behavior of the membranes is also inves- 
tigated analytically. 

1 Introduction 

Membranes serve a number of general functions in our cells and tis- 
sues. They separate cells and cell compartments [I112]- They act as 
barriers as their lipid core is permeable to water and to small hydro- 
carbon molecules and relatively impermeable to macromolecules and 
polar molecules. In addition, membrane proteins facilitate or assist 
the transport of ions and macromolecules into and out of the cells. 
Since biomembranes play a vital role in biological process, they have 
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attracted significant experimental and theoretical interest [3l [H [5l [6] . 
The physics of membrane adhesion is one of the aspects that has got 
considerable attention [7l[8l[9l[l0l[IIl[l2l[l3l[Ill[l5l[l6l[I7]. Re- 
cent experimental realization [18] uncovers the formation of domains 
between short and long receptors of T cells. The dynamics of adhe- 
sion induced phase separation was also considered theoretically in the 
works |191I20| and its detailed equilibrium studies were reported in the 
previous studies |22ll23t [2^. All these theoretical and experimental in- 
vestigations have significantly contributed to the basic understanding 
of the physics of membranes. More recently, we presented a theoreti- 
cal study for adhesion-induced lateral phase separation of membranes 
with short stickers, long stickers and repellers confined between two 
hard walls [25|. The effects of confinement and repellers on lateral 
phase separation were investigated. It is found that the sticker criti- 
cal potential depth tends to increase as the distance between the hard 
walls decreases which suggests confinement-induced or force-induced 
mixing of stickers. 

Despite the success to understand the thermodynamic properties 
of these non- homogeneous membranes, more efforts are needed to get 
a complete picture of the physics of membranes due to the fact that 
membranes contain large number of macromolecules which are orga- 
nized in complex fashion. Most of the recent works considered mem- 
branes interacting with one or two types of stickers (receptor /ligand 
pairs). However the adhesive molecules may have different resting 
length; the experimental investigations of T cells unveiled two or more 
receptor-ligand pairs being involved in the binding process [15]. An- 
other crucial issue on the adhesion induced lateral phase separation is 
the effect of confinement. Cell adhesion often occurs in the presence of 
external force field due to external flow [25] . The aim of this work is to 
investigate the phase behavior of confined membranes interacting with 
more than one types of stickers. First we investigate the phase behav- 
ior of membranes with short and long stickers analytically which can 
be taken as an independent check of the results found in the previous 
works (Ml [25] . We then extend our previous studies by considering 
membranes interacting with three adhesive molecules. 

The rest of the work is structured as follows. In section II, we 
introduce our model. In section III, we present the mean field approx- 
imation while section IV and V deal with membranes interacting via 
two and three types of stickers, respectively. Section VI deals with 
summary and conclusion. 
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2 The model 



We consider non- homogenous membrane that interacts via receptor /hgands 
(stickers) of different characteristic length as shown in Fig. 1. For sim- 
phcity, a discretized membrane (lattice gas model [121 121] ) with lattice 
constant a is considered. The separation field / separates the mem- 
brane from the membrane (substrate) while the occupation number 
nj = shows no sticker present at the lattice site i while rii = 1, 2, 3... 
indicates the presence of stickers 1, stickers 2, stickers 3... at lattice 
site i, respectively (see Fig. 1). 




Figure 1: Schematics diagram of membrane with stickers of different size 
near to a substrate, z denotes the local separation field that separates the 
membrane from the substrate. 



Introducing the rescaled fields z = {l/a)^Jk/T and V = V/T re- 
duces the number of independent parameters. For tensionless non- 
homogenous multi-component membranes, the Hamiltonian of the 
model is given by 

H[z,n] = H^i[z\+^5i,rH{yi{zt)- + 
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with the elastic term 
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(1) 
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where H^i [z] denotes the discretized bending energy of the membranes 
with an effective bending rigidity k. Typically, k = 10 — lOksT. The 
discretized Laplacian is given by A^Zj = zn + Zi2 + Zis + -Zj4 — 4zj. 
Vi{zi), V2{zi) and V3{zi) denote the potentials for stickers 1, stickers 
2, and stickers 3, respectively while fii, ^2 and /is denote the chemical 
potentials for the stickers 1, stickers 2 and stickers 3, respectively. The 
partition function Z can be written as 



n 



dzi 



W exp[-i7[z,n]] 



i ni=0,l,2 

Tracing out the sticker degrees of freedom one gets 
Z 



(3) 



n 



dzi [eyi-p[-Hei{z) 



(1 + exp [-Vi{zi) + i^Li] + exp [-V2{zi) + /is] + 
exp [-V^{zi) + /i3] + ...)• 

Rearranging terms, Equation (4) can be written as 

Z = (1 + exp[/ii] + exp[/Z2] + exp [/is] + exp [/is] 



(4) 



n 



dzi exp 



-H,i{z) + Y,[V^ff{z,)] 



(5) 



The effective potential V^f^{z) is given by 

V^ff{z) = -ln[(l + exp[-Vi(z) + /zi] + exp[-y2(^) 
exp[-y3(2)+M3] + -)/Co] 



M2] + 



(6) 



where ^0 = 1 + exp[/ii] + exp[/i2] + exp[/i3] + .... Integrating out the 
stickers degrees of freedom leads to a homogeneous membrane with an 
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effective potential which is given by Eq. (6). The stickers potential is 
modeled as square- well potential [121 121]) Vi = Ui for zi < Zi < Z2, 
V2 = U2 for zz < Zi < Z4, V3 = for z^ < Zi < zq and Vi, V2, 
V3...=0 otherwise. The effective potential (6) then can be rewritten 

as V^ff = 00 for z < 0, V^^^ = U^^^ for zi < Zi < Z2, V^^f = U^^^ 
for z,<z.< Z4, V^ff = Ul^f for z, < z. < z, ... where 

Ul^^ = - In [(1 + exp[-C/i + /ii] + exp[//2] + exp[/i3] + ..O/^o], (7) 

U^f^ = - In [(1 + exp[-C/2 + /i2] + exp[/ii] + exp[/X3] + ..O/^o] (8) 
and 

U^ff = _ln [(1 + exp[-[/3 + ^3] + exp[^i] + exp[/i2] + (9) 

One should note that 23 > 2:2, 25 > Z4 and so on, i.e. the binding wells 
do not overlap. 

In the next section, we study the phase behavior of membranes 
with short and long stickers utilizing the mean field approximation. 



3 Mean-field approximation 

When fluctuation of the membranes not too strong, one can apply 
mean-field approximation to the discretized Laplacian A^Zj = 4[zmm — 
Zi\ where Zmin designates the average separation field. Substituting 
this equation in Eq. (1), the elastic term takes a simple form: H[z\ = 
'Ylii{'^[^min — Zi])^ . Within the mean field approximation, one can write 
Eq. (5) as 



Z 



(Co) 



N 



dzi exp I 



(10) 



The free energy of the membrane after some algebra is given by 



G 



ln[ 



dzi [exp [ 



/eo]. (11) 



Varying the free energy (11) with respect to Zmin leads to a self con- 
sistence equation of the form: 



ze^^[-S{zmin - Zif - yf^iiyidzi 



rOO r 

Jo exp[ 



Zif - V-}f{i)]dz, ■ 



:i2) 
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Figure 2: Double-well potential V^^^ with two degenerate minima at both 
sides of the wells. The potential has two square wells of depth Ui'^^ and 
U2^^^ within the range Zyj^. = Z2 — Zi and z^^ = Z4 — 23, respectively. The 
effective potential V^^-^" = 00 for 2 < and z > z^. We fix Za = z^ — z^ = Zi. 



3.1 Membranes with short and long stickers 

Let us now consider membranes interacting via two types of stickers 
whose equilibrium phase behavior are governed by an effective double- 
weh potential with the effective depth U^^-^ and U2^-^ as shown in Fig. 
2. U^^^ and U2^-^ depend on the stickers binding energies Ui and 
U2, respectively. Hereafter, all energetic quantities are given in unit 
of ksT and to make the model analytically solvable (for symmetry 
reason) let zi = 0, Z4 = 2:5, Zwe = Z2 — zi = — Z3 and Zba = Z3 — Z2- 
When membranes are confined in one of the wells of the effec- 
tive double-well potential, they remain stable with respect to thermal 
fluctuations provided that the potential wells are deep enough, mem- 
branes exhibit two coexisting states with two average separations Zmin 
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Figure 3: Graphical solution R = zexp[—8{zmin — — 

V^^^ {i)]dzi/ Jq°° exp[— 8(2;m,m — ZiY — V^^^ {i)\dzi versus Zmm for parameter 
choice Zwe = 0.4 and Zba = 0.4. The dashed and the thick solid lines repre- 
sent the graphical solution for [[/^■^■^l = 8 and IW^-^-^ = 0.2, respectively. 



[3] . In this case the membranes remain trapped in one of the wells and 
exhibit a first order transition. However for weak potential wells, the 
membrane can easily surmount the potential barrier and its average 
separation field Zmin = 0. 

At what critical potential depth Uq^"^ do membranes confined in 
the double- well potential "tunneF? This can be addressed by solving 
Eq. (12) numerically. For instance for parameter choice Zy^e = 0.4, 
Zba = 0.4 and = 8, the graphical solution is depicted in Fig. 3. 

The figure clearly reveals Zmin attending two distinct values z^in = 0.2 
and Zmin = 0.98 which exhibits llf^-^^l > \Uq''^\. On other hand 
for Zwe = 0.4, Zba = 0.4 and llf^^^ = 0.2, z^m = revealing 
\U^ff\ < \Ul/^\ (see Fig. 3). Furthermore for parameter values of 
Zwe = 0.4 and Zba = 0.4, the dependence of z^m as a function of W^f^ 
is plotted in the Fig. 4. The same figure depicts when llf^-^-^] > \U^^'^\, 
membranes are confined within one of the potential wells. As the ef- 
fective potential depth decreases, the tendency for membranes to stay 
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in the potential wells decreases. When \U'^-I'^ < \U^r \, membranes 
overcome the barrier height and as the result Zmin = 0. 

The critical potential depth , below which phase separation 
occurs, can be obtained by minimizing the free energy, i.e.; —d'^G/dz^^^ = 
0. For fixed zi = and Z4 = 25, one finds the critical potential depth 

Ul^f = ln[l - eM-KZu,e{Zu,e + ^ba)](l + ^)]. (13) 

Zba 

Exploiting Eq. (13), one confirms that when the potential width Zy^^ 
steps up, declines which is consistent with the work |24t l25|. Note 
that Zfea signifies the length difference between the two stickers. As 
it can be readily seen in Eq. (13), when increases, declines. 
This is because when the length difference between the two stickers 
increases, the energy cost of keeping the short and long stickers in 
close proximity increases. 

The effective critical potential depth Uq^ relies on sticker binding 
energy U and stickers chemical potential ^. Comparing Eqs. (6) with 
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(13) and assuming = fJ.2 = one obtains the stickers critical 
binding energy 

Uc = H-H^] (14) 

where 

= eM+8.»e(.»e+.fc.)^^^ + (l + gM)(^^^ + 2z^e) (15) 

and 

/?2 = -l + e^'-^^'-^+"^-hba-2z^e. (16) 

When the stickers binding energy Uc is less than the critical binding 
energy U < Uc, the short and long stickers segregate into distinct 
domains while when U > Uc, the stickers become mixed. The depen- 
dence of Uc on model parameters can be explored via Eq. (14). When 
the length between short and long stickers Zba as well as the width of 
the stickers potential increases, the two types of stickers separate 
into two distinct domains at smaller critical binding energy Uc- 

Let us now consider a case where zi ^ 0, zi = z^ — z^ = Za (the 
two wells of the effective double-well potential are far from the hard 
wall) as shown in Fig. 2. In this case the effective critical potential 
depth Uc'^'^ after some algebra is given by 

^/^•^■^ = ln[exp[(8Za?n'l)]((exp {8Zy,e{Zwe + Zba)) - 'i^)zba - 

2z^e)/a] (17) 

where 

a = Zba + Zba exp[(8(2;a + Zu,e)m2] - exp[(8zami)]zfea + 

2{Za + Zi,e- ex.p[{8za'mi)]zwe)- (18) 

The values for mi and m2 are given by ml = Zba + Za + 2zwe and 
m2 = Zba + Za + Zwe- One can note that in the limit Za — ?• 0, Eq. 
(17) converges to Eq. (13). Utilizing Eq. (17), one can see that as 
Zwe and Zba increase, Uc^-^ declines which suggests that as the width 
of the stickers potential as well as the length between the short and 
long stickers increases, the stickers segregate into two distinct domains 
at a lower stickers critical potential depth. In the limit Zba oo 
or Zwe oo, \Uc\ while in the limit Zba — )> or z^e 0, 
\Uc\ oo. When Za increases, the entropic repulsion of membranes 
with the wall reduces. Thus the critical potential depth of the stickers 
Uc"^"^ decreases which exhibits that confinement facilitates demixing 
of stickers and hinders lateral phase separation. 
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4 Membrane with short, intermediate 
and long stickers 

Let us now consider membranes with short, intermediate and long 
stickers near to the substrate. Eliminating the stickers degree of free- 
dom results in membranes interacting with triple-well potential with 
the effective depth V^^^ , ^2^^ and Ut^^^ as shown in Fig. 5. U^^^ ■, 
U^^^ and Ul^^ are function of the stickers binding energies Ui , U2 and 
^73, respectively. Once analyzing the effective critical potential depth, 
one can retrieve the critical point of the stickers binding energy. 
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Figure 5: Triple-well potential V^-'-' with three degenerate minima. The 
potential has three square wells of depth Ui^^^, t/a"^^ and U-s^^^ within the 
range = 2:2 — zi = Z4 — = zq — z^, respectively. Here Za = zi = zj — zq. 
The effective potential V'^^-l' = 00 for z < and z > zj. 



Our earlier analysis unveils that for membranes with two types 
of stickers, the length disparity between the short and long stickers 
causes lateral phase separation. For pronounced length difference. 
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lateral phase separation occurs even at high temperature or equiva- 
lently at shallow critical potential depth. However the situation is 
quiet different when there are intermediate stickers embedded in the 
membranes. The presence of intermediate stickers reduces the bending 
energy cost as the length gap between the short and long stickers gets 
diminished. Hence in the presence of additional intermediate stickers, 
the lateral phase separation among the stickers occurs at lower tem- 
perature or deeper critical potential depth which suggests that lateral 
phase separation occurs at the expense of high sticker binding energy 
and thus the intermediate stickers endorse mixing of stickers. 




Figure 6: The critical effective potential depth U^[; = U^J versus t/^' 
parameter choice = 0.1, = = 0.16. 

The dependence of the critical potential depth on the model pa- 
rameters can be explored numerically via Eq. (12). In order to con- 
ceive the effect of the intermediate stickers, let us vary 1/2^ and for 
simplicity let U^^-^ = U^^^ . Figure 6 depicts the plot U^j^^ = U^l^ as 
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a function of C/2 parameter values — 0.1, z^a — -^we 

= 0.16. 

When U^^^ = 0, \U^c^\, \U^[^ = 1-326 while as {U^^^ steps up, 
l^ic'^l; \^3c^\ increases revealing the intermediate stickers hindering 
the lateral phase separation between the short and long stickers. Fig- 
ure 7 exhibits the plot of |f/if'^|, \U^c'^\ as a function of z^a for fixed 
Ui'^'^ , Za = 0.1 and z^e = 0.16. As zi,a increases \Ull^\, \U^l^\ de- 
clines; when \U2^^\ steps up iC/^f'^l, increases once again shows 
that the intermediate stickers facilitate mixing of the stickers. 

These ideas can easily checked in experiment by considering vesicle 
adhesion to supported membrane with several types of stickers. Our 
theoretical prediction reveals that the presence of intermediate stick- 
ers suppresses lateral phase separation. This effect can be checked 
by incorporating different types of adhesive molecule of different size 
to membrane surface. As the number of species of stickers increases, 
lateral phase separation gets suppressed. However lateral phase sep- 
aration of two or more types of stickers has been observed in experi- 
ment. This leads to the conclusion that either nature keeps its types 
of stickers to a limited number or equilibrium modeling, although it is 
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simple, fails to capture the phase behavior of membrane with several 
types of stickers. 

5 Summary and conclusion 

The equilibrium phase behavior of membranes with stickers of several 
types is explored at the mean field level. We first consider membranes 
with short and long stickers. We investigate how the critical potential 
depth \Uc'^-^\ behaves as the function of the model parameters. The 
central result shows that \Uc'^^\ decreases not only as the width and 
length of the stickers increase but also when Za increases. In the 
presence of intermediate stickers, the numerical results uncover that 
the intermediate stickers hinder the lateral phase separation. 

In conclusion, via mean field approximation, we study adhesion in- 
duced lateral phase separation of membranes interacting with multiple 
species stickers of different sizes. Though the mean-field approxima- 
tion neglects the effect of fluctuations, we believe that the qualitative 
behavior of the result obtained will not be affected and thus this work 
is crucial not only for fundamental understanding of the physics of 
membranes but also for the construction of artificial membranes. 
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